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In this Letter, the dynamic phase transitions of the time-dependent Ginzburg-Landau equations 
are analyzed using a newly developed dynamic transition theory and a new classification scheme of 
dynamics phase transitions. First, we demonstrate that there are two type of dynamic transitions, 
jump and continuous, dictated by the sign of a nondimensional parameter R. This parameter is 
computable, and depends on the material property, the applied field, and the geometry of domain 
that the sample occupies. Second, using the parameter R, precise analytical formulas for critical 
domain size, and for critical magnetic fields are derived. 

PACS numbers: 74.20.-z, 74.20.De, 74.25.Dw 



One central problem in the theory of superconductiv- 
ity is the dynamical nature of the phase transition be- 
tween a normal state, characterized by a complex or- 
der parameter ip that vanishes identically, and a super- 
conducting state, characterized by the order parameter 
that is not identically zero. In this Letter, we address 
this problem by conducting rigorous theoretical analy- 
sis on dynamic phase transitions for the time dependent 
Ginzburg-Landau (TDGL) model. 

TDGL model. Let tt C K" (n = 2, 3) be a bounded 
domain, ip : SI — > C the order parameter, H : Vl — > 1" the 
magnetic field with A being the magnetic potential given 
by H — curl A, and H a the applied field with potential 
A a such that curM a ~ H a . The nondimensional TDGL 
equations for (ip, A, </>), with A being the deviation from 
the applied field, are given by: 
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Here the nondimensional parameters are defined by 

a = -2aVbm s D/e 3 s h, (3 = 2m s D/h, 

fi^hD/Vbe s , ( = 4no-le 2 /c 2 h 7 

7 = ^ire 2 /m s c 2 l, A = A(T) = {m s c 2 b/ATre 2 s \a\) 1/2 , 

k = A/£, i] = 4naD/c 2 , 



where h is the Planck constant, e s and m s are the charge 
and mass of a Cooper pair, c is the speed of light, 
A = A(T) is the penetration depth, £(T) is the coher- 
ence length, r for the relaxation time, and k = A/£ is the 
Ginzburg-Landau parameter. 

The TDGL equations are supplemented with an initial 
condition for (ip, A), a free-slip boundary condition for 
A, and either the Neumann or the Dirichlet or the Robin 
boundary conditions for i/j; see de Gennes Here We 
shall see in later discussions that the parameter a plays 
a key role in the phase transition of superconductivit 
which is given in terms of dimensional quantities by 



i 



= a(T) = 



2Vbm s DN T c -T 



(2) 



where TVo the density of states at the Femi level, and T c 
the critical temperature where incipient superconductiv- 
ity property can be observed. 

Dynamic transition. As mentioned earlier, the 
study of dynamic phase transition problem for the TDGL 
equations is based on a new dynamic transition theory 
by the authors @, 0] . The starting point of the theory 
is to put the TDGL equations into the perspective of an 
infinite-dimensional dynamical systems as follows: 



— = L\u + G(u, A), 
at 



u(0) = uq. 



(3) 
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where u : [0, oo) —tHis the unknown function, H is a 
Hilbert space, L\ is a linear operator, G(it, A) is a nonlin- 
ear operator, and A is the system parameter. Then under 
proper physical conditions, the dynamical transitions of 
a basic state of © at a critical parameter Ao, where 
the principle of exchange of stabilities hold true, can be 
classified into three categories: continuous as shown in 
Figure [TJ jump as shown in Figure [U and mixed; see 
the above references for further details. This theory has 
been applied to the TDGL equations ([T]) to characterize 
the dynamic transition, will be be used in the analysis 
hereafter in this Letter. 
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Eigenvalue problem: Consider the following equa- 
tion 



(i/jV + A a ) 2 ip = atp Vx £ fi, 



(4) 



with one of the boundary condition for ip. There are an 
infinite real eigenvalue sequence of (|4]) as < a\ < a 2 < 
■ ■ ■ < ecu < • • • — * +oo. The eigenvalues of ((4|) always have 
even multiplicity. In this Letter, we consider only the 
case where the first eigenvalue ot\ of ^ has multiplicity 
two, i.e., ot\ is complex simple eigenvalue. We let 



e = e± — ipn + itp 



12. 



C2 



-ipu + iip 



ii- 



(5) 



be the eigenvectors corresponding to ol\. 

A necessary condition: In superconductivity, the 
parameter a in P can not exceed a maximal value 
a(T) < a(0) because T > 0. Hence, a necessary con- 
dition for the possible phase transition from the normal 
state to superconducting states is: 



at < a(0) 



2Vbm s DN 



(6) 



where a\ is the first eigenvalue of ([4]). 

For the case where the Neumann boundary condition 
dip/dn = on dil, the first eigenvalue a\ — 0, which is 
independent of fl, the geometry of the sample. Therefore, 
the condition ^ always holds true. 

However for the Dirichlet and the Robin boundary con- 
ditions, the situations are different. We know that the 
first eigenvalue a,\ depends on Q. We have for example 
lim|Q|^o a i — °°- Hence, the condition ([6]) implies that 
for the cases with either the Dirichlet or Robin bound- 
ary conditions, including the case where the sample is 
enclosed by a magnetic material or a normal metal, the 
volume of a sample must be greater than some critical 
value |fi| > V c > 0. Otherwise no superconducting state 
occurs at any temperature. 

Physical theory and experiments show that there is 
a critical applied magnetic field H c by which a super- 
conducting state will be destroyed, and H c satisfies the 
following approximate equation near the critical temper- 
ature T c which is given by © 



H C (T) = H Q (l-T 2 /T 2 ] 



(7) 



Equation ([7]) is an empirical formula. A related equa- 
tion is the equation of critical parameter: 



a(T) = ai (H a ,Q), 



(8) 



where a(T) is the parameter in ([T]) and a\ = ai(H a ,Q) 
the first eigenvalue of Q . It is expected that the applied 
magnetic field H a satisfying ((8|) is the critical field H c . 

Definition of R: We start with the introduction of a 
crucial physical parameter, which completely determines 
the dynamical properties of the phase transition of the 
Ginzburg-Landau equations. 



For the first eigenfunction e, the following Stokes prob- 
lem has a unique solution: 



curl 2 _4 + = \e\ 2 A a + ^i(e*Ve - eVe*), 
divA = 0, 

A • n\dn = 0, curL4 x n\ dn = 0. 
Then we define a physical parameter R as follows 



(9) 



R = - 



/3 2/ n |curL4o| 2 <fe 



dx 



(10) 



From ([9]) and (fT0|) it is easy to see that the parameter 
R is independent of the choice of the first eigenvectors of 
dH) and H — curL4 given by © depend on the applied 
magnetic potential A a and the geometric properties of 
f2, the parameter R is essentially a function of A a , and 
physical parameters j3, 7, /i. 

As mentioned earlier, there are two phase transitions: 
Type-I and Type-II, determined by a simple parameter 
R defined by (|T0|) . This parameter R links the supercon- 
ducting behavior with the geometry of the material, the 
applied field and the physical parameter's. 

An equilibrium state (tp, A) of the TDGL equations ^ 
is called in the normal state if ip = 0, and it is called in 
the superconducting state if ip ^ 0. A solution (ip,A) of 
(fTj) is said in the normal state if (ip,A) is in a domain of 
attraction of a normal equilibrium state, otherwise (ip, A) 
is said in the superconducting state. 

Dynamic phase transition for R < 0: By the phase 
transition theorems obtained by the authors [J] , the crit- 
ical temperature of superconducting transition satis- 
fies T\ < T c where T c is as in © and T = T 1 satisfies ®. 
It is known that lim| J 4 a |^ 00 ai(H a , £1 = 00, which implies 
that the applied magnetic field H a can not be very strong 
for superconductivity as required by condition ([B]). 

We have shown in [3| that when a > oti, the equa- 
tions |T]) bifurcate from ((tp, A), a) = (0, ai) to a cycle of 
steady state solutions (ip a iA a ) which is an attractor at- 
tracting an open set U\T C H. From the physical point 
of view, this theorem leads to the following properties of 
superconducting transitions in the case where R < 0: 

First, when the control temperature decreases (resp. 
increases) and crosses the critical temperature T 1 , there 
will be a phase transition of the sample from the normal 
to superconducting states (resp. from the superconduct- 
ing to normal state). 

Second, when the control temperature T > T 1 , un- 
der a fluctuation deviating the normal state, the sample 
will be soon restored to the normal state. In addition, 
when T < T^, under a fluctuation deviating both the 
normal and superconducting states, the sample will be 
soon restored to the superconducting states. 

Third, in general, the supercurrent given by 



J s (a) = — j(A a + A a )\ip a f 



i(lp* a Vlp a - ifraVlpa) 



is nonzero, i.e., J s ^ for T < T* (a\ < a) 
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Fourth, the order parameter ip a and supercurrent 
J s (a) depend continuously on the control temperature 
T (or the parameter a): 

4> a ^0, J»->0 ifT-» (T c x )- (or a^a+). 

Fifth, the superconducting state of the system is dom- 
inated by the lowest energy eigenfunction of in the 
sense that ip — C\a — a\\ l l 2 e + o(\a — ail 1 / 2 ). 

Sixth, the phase transition is of second order in the 
Ehrcnfest sense with the critical exponent (3 = 1/2, and 
its phase diagram is as shown in Figure [TJ 




FIG. 1: The continuous transition with R < 0. 

Dynamic phase transitions with R > 0: Consider 
a material described by the TDGL model with R > 0. 
There are two transition temperatures T® and (Tj? > 
T}) such that 

ot(Tl) = cti (i — 0, 1) with ao < ai, 

and the following phase transition properties hold true: 
First, when the control temperature T decreases and 
crosses or equivalently a increases and crosses a±, 
the stability of the normal state changes from stable to 
unstable. 

Second, when T* < T < T° (a < a < ai), physi- 
cally observable states consist of the normal state and the 
superconducting states in E 2 (see Figure 8.20). When 
T < {ot\ < a), physically observable states are in E 2 . 

Third, when the control temperature T is in the inter- 
val T]. < T < T c ° (or «o < a < Oil), the superconducting 
states in E" are unstable, i.e., with a fluctuation devi- 
ating a superconducting state in E", transition to either 
the normal state or a super-conducting state in SJ? will 
occur. 

Fourth, at the critical temperature T c ° (resp. at T*) 
of the phase transitions, there is a jump from the super- 
conducting states to the normal state (resp. from the 
normal state to superconducting states). 

Fifth, the other energy-level eigenfunctions possibly 
have a stronger influence for the superconducting states. 

Sixth, in the temperature interval T c x < T < T c °, 
phase transitions occur and are accompanied with the 
latent heat to appear. 

Based on the conclusions (l)-(4) above, we can draw 
the phase diagram in Figure [21 where the critical temper- 
ature T c in the interval T\ < T c < Tj? is the transforming 
point. 




FIG. 2: The jump transition with R > 0. 



Critical Sample Size and Critical Magnetic 
Fields: Theories and experiments illustrate that the ge- 
ometry of samples 51 and applied magnetic fields H a have 
important influences for the superconducting behaviors. 
In the following, we shall apply the formula (|10|) and 
eigenvalue equation (01 to discuss this problem. 

For < L < oo and h > 0, let Q = D x (0, h) C K 3 , 
x' = (xi,x 2 ) € -Do, G Do, and 

Q(L) = {(Lx',x 3 )\ x' eD a ,0<x 3 < h}. 

Let the applied field be given by H a = H a (x') = 
(0, 0, H(x')), where x' — (xi,X2) £ Dq. Then H a induces 
an applied field H a on Cl(L) by H(y) — H{y/L) for any 
y = Lx' with x' £ Dq. 

Let H a = curL4 a , i.e., H = ^ - §§f , where 

A a = (Ai(x'),A 2 (x%0). (11) 

Then we can show that the nondimensional parameter 
R = R(L, H) on Q(L) is given by 

R{L,H) = - K 2 fi 2 + 2L{p 3 L i +p 2 L 2 + Pl ), (12) 

where for any < L < oo, 

( f n |curL4 |W 

J Do \ e \ dx 
K curl 2 A + V0 = \e\ 2 A a in x £ D , 
( _ V^JcurL9op^ 

[ curl 2 Do + = ei2Ven in x' £ Do, 
4[i j Do curlAo • curlDo<ia;' 
P2= ° L kl 4 ^' ' 

Here 5 > is independent of L, and e(x') = en + iei2 is 
the first eigenfuction of 

(i/uV + L 2 A a ) 2 e = L 2 ae in D , 

de (13) 

— = on dD . 

on 

Formula (|12p has many applications. In particular, we 
derive the following effort of the domain and the applied 
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field. It is easy to observe that p% is essentially a |if| 2 
term, p 2 is a |-ff |, and p± is a zeroth order of \H\. 

First, by ([T2|l . we derive that given a superconducting 
material, and an applied field H a =/= 0, there is a critical 
scale Lo > such that the phase transition in fi(L) is a 
continuous transition if L < Lo, and is a jump transition 
if L > Lq. In addition, Lq is the unique real root of 

(XsLq + a 2 Ll + a\L = -n' > 5 '. 

This theoretical conclusion is known experimentally for 
the small sample Vt(L) (i.e., 1 3> L) case. 

Second, for Type I superconducting materials (k 2 < 
1/2), there are three critical magnetic fields H Cl < H C2 < 
H C3 such that the following hold true: 

1. If < H < H Cl , then the transition is Type-I and 
the superconducting state is in the Meissner state; 

2. If H Cl < H < H C2 , then transition is Type-II and 
in the Meissner state, 

3. If ai(A C2 ) = a(0), then for any H C2 < H, it is in 
the normal state. 

4. If ai(A C2 ) < a(0), then 

(a) if H C2 < H < H C3 , then transition is Type-II 
and in the mixed state, and 

(b) if H C3 < H, then it is in the normal state and 
ai{A C3 ) = a(0). 

Third, for type II superconductors, the results are sim- 
ilar: 

1. If < H < H Cl , then the transition is Type-I and 
in the Meissner state, 



2. If H C1 < H < H C2 , then transition is Type-I and in 
the mixed state, 

3. If a\{A C2 ) = a(0), then for any H C2 < H, it is in 
the normal state. 

4. If ai{A C2 ) < a(0), then 

(a) if H C2 < H < H C3 , then transition is Typc-II 
and in the mixed state, and 

(b) if H C3 < H, then it is in the normal state and 
ai(A C3 ) = a(0). 

We remark that the properties above are well known 
by the Abrikosov theory and physical experiments. How- 
ever, here conclusions are derived rigorously using the 
recently developed dynamic transition theory, and the 
critical fields can be solved from (fl"2"|) . 

Fourth, for a small sample, there is a remarkable dif- 
ference that the first eigenvalue a\{A a ) of with the Neu- 
mann boundary condition is bigger for small applied 
magnetic field H a = curL4 a . Hence the critical mag- 
netic field H c = curL4 c satisfying that a\{A c ) = a(0) is 
very small, which implies that for a small sample, the 
superconductivity is only in the Meissner state. 

Finally we would like to mention that there have been 
extensive studies on bifurcation and stability analysis for 
superconductivity; see among others [alsHjIS 9]. The 
study in this article is based on a newly developed dy- 
namic transition theory by the authors 0, Q . With this 
new theory, many long standing problems in phase tran- 
sition problems in science and engineering are becoming 
more accessible. In particular, applications are made for 
various models from science and engineering, including, 
in particular, problems in statistical physics, classical and 
geophysical fluid dynamics. 
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